Energy transfer across magnetopause under dawn–dusk IMFs

A parametric study on the energy transfer of the solar wind across the magnetopause entering the magnetosphere is conducted using a global magnetohydrodynamic numerical simulation. The characteristics of the mechanical and electromagnetic energy distribution under the dawn–dusk interplanetary magnetic fields (IMFs) are investigated by analyzing magnetic reconnection and viscous effect, and compared with the radial and north–south IMFs. It is shown that (1) the interactions at the magnetopause and the transfer of energy across this boundary move in relation to the IMF orientation. (2) For the duskward IMF, the mechanical energy flow clearly enters the equatorial and low-latitude regions on the dayside, and the electromagnetic energy flow has a small inflow on the equatorial and low latitudes of the dayside. A significant energy inflow appears on the dawn side in the northern hemisphere and the dusk side in the southern hemisphere near the polar cusp. (3) The energy distribution characteristics across the magnetopause under dawn–dusk IMFs are mirror symmetric about the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y=0$$\end{document}Y=0 plane. (4) For a magnetic field of 5 nT, the electromagnetic energy input under the dawn–dusk IMFs is twice as large as the mechanical energy and the electromagnetic energy under the radial IMF, which is five times as large as the electromagnetic energy during the pure northward IMF, but only half as large as the electromagnetic energy under the pure southward IMF. The mechanical energy input under dawn–dusk IMFs has the same magnitude as that under radial and north–south IMFs. The magnitude of the energy transfer rate for the dawn IMF and dusk IMF (about 3.5%) is between 1.71% for the northward IMF and 4.95% for the southward IMF, but higher than 2.22% for the radial IMF. The Akasofu-type energy-coupling formula, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon$$\end{document}ε, underestimates the energy input from the solar wind under \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_{y}$$\end{document}By dominated IMF.

The magnetic reconnection between the interplanetary magnetic field (IMF) and geomagnetic field is the most critical process that controls the energy transfer from the solar wind to the magnetosphere across the magnetopause. Compared as the geomagnetic field is relatively stable, the variable IMF as a part involved in the interaction naturally plays a key role in this coupling system. The reconnection can occur for nearly any orientation of the IMF and produce different reconnection rates under different IMF directions. The reconnection efficiency affects the amount of energy input 1 , therefore, understanding the energy distribution characteristics, the amount of energy input, and energy transfer efficiency under different directions of IMF is worth exploring.
The effect of the north-south components of IMF, B z , on the magnetopause magnetic reconnection of the magnetosphere is the most pronounced, and many works have investigated the magnetic reconnection and the energy transport under the northward or southward IMF [2][3][4] . It was demonstrated that the reconnection for southward IMF always occurs near the equatorial plane on the magnetopause, where the reconnection rate is the greatest. In contrast, if the IMF is northward, the magnetic reconnection always occurs near the tail behind the polar cusp, and the reconnection becomes the weakest 5 . Using the polar cap potential difference as an equivalent proxy for the magnetic reconnection rate, the reconnection rate under northward IMF is about 10-30% of that under southward IMF 6 .
The impact of the radial component of the IMF on magnetopause 7 , magnetic reconnection, and energy transfer across the magnetopause was investigated recently. A radial IMF can make the reconnection location and energy transfer asymmetric in the northern and southern hemispheres 8,9 . Petrinec et al. 10 found that the reconnection site on the magnetopause during radial IMF (without considering the strong radial IMF) is not exactly at the theoretically estimated position but is located further downstream of the magnetopause. Using a global magnetosphere MHD model combined with the THEMIS observations, Pi et al. 11 studied the global structure of the adjoining regions near the magnetopause under radial IMF and showed that there is a hemispheric asymmetry of the reconnection location of the dayside magnetopause. For a non-null Earthward IMF B x , the magnetic reconnection occurs on the dayside magnetopause in the southern hemisphere and behind the cusp in the northern hemisphere. Similarly, using a global MHD simulation, Peng et al. 12 found that the dayside reconnection line of the magnetopause moves northward (southward) when the B x 0 ( B x 0) , resulting in the magnetopause movement. Hoilijoki et al. 13 used another global MHD simulation to find that the radial IMF mainly affects the reconnection position on the dayside magnetopause and has a more significant impact on the intensity of energy transfer in the load zone (the region where energy is converted from magnetic to kinetic form, such as the subsolar point of the magnetopause where reconnection occurs.) Lu et al. 14 used a global MHD model to study the energy transfer across the magnetopause under radial IMF and found that the energy transfer across the magnetopause is asymmetric and can not be estimated by the Akasofu-type energy-coupling formula. The formula is ε = 4π where v is the solar wind velocity, B is the IMF intensity, L 0 = 7 R E is an empirical scaling parameter, and θ is the IMF clock angle ( tan θ = B y B z ). Particularly, the electromagnetic energy and the energy transfer rate of radial IMFs is lower than the one of pure southward IMF but higher than the one of pure northward IMF.
The dawn-dusk component of the IMF also affects magnetopause and magnetic reconnection 15,16 . When the Y component of the IMF dominates, the interplanetary magnetic field is reconnected with the geomagnetic field on the day side. The open magnetic field lines connected to the Earth move to the magnetotail and reconnect with the geomagnetic field near the tail lobe of the magnetosphere again, forming the low latitude boundary layer 17,18 . Marcucci et al. 19 found that the direction of the accelerated plasma flow changes when the IMF turns from dusk to dawn, as evidence of a change in the reconnection position at the magnetopause. Park et al. 20 calculated the electric field at the magnetopause and found that anti-parallel reconnection is more important than the component reconnection for the case of slight dipole tilt and IMF B y component. Komar et al. 21 tested and compared a number of existing models predicting the location of magnetic reconnection at dayside magnetopause for various solar wind conditions. They found that the maximum magnetic shear model does not rotate with the magnetic separators for different IMF clock angles in the simulation without dipole tilt. Cowley et al. 22 discussed the observed magnetospheric asymmetry caused by IMF B y according to an open model of the magnetosphere. In addition, the dawn-dusk IMFs can also lead to the asymmetric distribution of the field aligned current (FACs) in the reconnection region of the magnetopause and increase the FACs in the reconnection region 23 . Using the Cluster data, Cheng et al. 24 selected 748 FACs cases and proposed a significant north-south asymmetry in the polarity of the FACs under the dawn-dusk IMFs, especially in the duskward IMF. A north-south asymmetry induced by the asymmetry of the tail lobe pressure occurs in the magnetosphere-ionosphere system when there is a significant IMF B y component, and the asymmetry appears during low tail reconnection and decreases during enhanced tail reconnection 25 . When the Earth's dipole is tilted in the direction corresponding to northern winter, the area of the polar cap under duskward IMF caused by the difference of the day-side reconnection rate is larger than that of dawnward IMF 26 . Cao et al. 27 proposed that the penetration with IMF B y into the magnetosphere is enhanced during strong magnetospheric convection. Guo et al. 28 used a global MHD simulation to investigate the effects of IMF B y on the FAC closure in the magnetosphere. The symmetry of the FACs across the noon-night meridional plane breaks in the ionosphere for dawn-dusk IMFs, unlike those cases in which the IMF is purely northward or southward. Using satellite observations and a global MHD model coupled with a ring current model, Holappa et al. 29 found that in the northern hemisphere, the flux of energetic magnetospheric protons and the growth rate of the circulation are larger for IMF B y 0 than for IMF B y 0. The above studies show that the magnetospheric structure and magnetic reconnection are affected by the IMF B y component, and the energy transfer characteristics across the magnetopause are also impacted. However, few quantitative studies on the distribution characteristics of energy transmission across the magnetopause and the energy transport rate have been reported under the dawn-dusk IMFs.
Another important factor contributing to mass and momentum transport across the magnetopause is the presence of viscous interactions between the solar wind and the magnetosphere that produce the magnetospheric convection 30 . It has been reported that the low latitude magnetopause is more conducive to the generation of K-H instability and viscous interaction when the IMF is northward 31 and the magnetotail viscosity under radial IMF is weaker than that in northward IMF 14 . However, the effect of viscous interactions under IMF B y is still unclear.
In this paper, through the global MHD simulations, we split the energy into mechanical energy and electromagnetic energy, and discuss their distribution on the magnetopause under B y dominated IMF, respectively. The energy transfer mechanism under dawn-dusk IMFs is also given by analyzing the location of magnetic reconnection occurrence and viscous effect. In addition, the energy transport rates under IMF B y are compared with the one obtained for the radial and north-south IMFs.

Global MHD model and methods
We use the Space Weather Modeling Framework (SWMF), a well-established set of computational models developed by the University of Michigan, to simulate the physical processes between the Sun and the Earth. SWMF is widely used to study the various effects of solar wind on the magnetosphere, such as the convection in the northward IMF 32-34 , the IMF B y 35,36 , dipole tilt angle 37 , and Parker-spiral IMF conditions 38 on the magnetosphere. The SWMF numerical model, including the fully coupled RCM module, is also used to study the influence of the southward component of the IMF 39 , the dynamics of magnetic storms 40 , and the energy transfer across the magnetopause 2,41 . The methods of our research are similar to the one we used to study the energy transmission across the magnetopause for purely southward and purely northward IMF 2, 14 .
In a typical simulation used in this work, the computational domain is defined by − 224 www.nature.com/scientificreports/ 64 R E , the grid size is 2 R E ; inside − 48 R E X 32 R E and − 32 R E Y, Z 32 R E , the grid size is 1 R E ; inside − 32 R E X 24 R E and − 16 R E Y, Z 16 R E , the grid size is 1 2 R E ; inside − 28 R E X 16 R E and − 12 R E Y, Z 12 R E , the grid size is 1 4 R E ; and inside − 14 R E X 10 R E and − 8 R E Y, Z 8 R E , the grid size is 1 8 R E . The inner boundary is a sphere with a radius of 3 R E .
To compare and discuss the influence of each of the IMF components on the energy transmission, we first analyze the influence of the variable IMF B y on the energy transmission across the magnetopause, a case that has not been studied before. As shown in Table 1, the solar wind conditions used in our study in cases 1-7 only change the strength and direction of IMF B y . In cases 8-11, we only changed the direction of IMF B x and B z , with a magnetic field strength of 5nT. The dipole tilt of all conditions is set to zero, and other physical quantities remain consistent.
This work investigates energy transfer under a steady solar-wind condition and IMF. As Kelvin-Helmholtz instabilities are not considered here, magnetic reconnection is then the major contributor to the coupling. Since we are dealing with the large-scale coupling of solar wind and magnetosphere, we did not make special set with resistivity.
For each case, SWMF is run for two hours to get a stable magnetosphere. Then, we adopt the method of tracking the solar wind streamline to automatically identify the smooth surface of the magnetopause and the approach of calculating the pressure gradient to determine the location of the polar cusp as done in previous studies 2, 14,34 . The intersection line between the posterior edge of the cusp (located at the maximum value of the pressure gradient) is located, and the identified magnetopause is defined as a boundary to separate the dayside and nightside magnetopause. Once the coordinates of the three-dimensional magnetopause surface are obtained, the transfer energy can be calculated 2,14,42 . The total energy input to the magnetosphere is defined as: where W m is the total mechanical energy input, W p is the total electromagnetic energy input, dS is the area of the surface element, and n is the unit vector perpendicular to the surface pointing outward the magnetopause. K m represents the mechanical energy flux density, which is oriented in the same direction as the magnetosheath flow field, and K p represents the electromagnetic energy flux whose direction is determined by the direction of E × B . The definition formulas of K m and K p are where P is the pressure, B is the magnetic field, v is magnetosheath velocity, E = B × v is the electric field, and the total energy density of the fluid element U is defined as U = P γ −1 + ρv 2 2 + B 2 2µ 0 ( γ = 5 3 ). Using the above equations, the mechanical energy and electromagnetic energy injected into the magnetosphere are calculated for all azimuths (the angle away from the positive Y-axis in the YZ plane) from the subsolar point to X = −40R E of the magnetotail in the YZ plane. When energy is injected into the magnetosphere, the energy value is negative, and vice versa.

Energy flow through magnetopause for dawn-dusk IMFs
Since the characteristics of energy transfer under the condition of B y = −5 nT and B y = 5 nT with P d = 3 nPa have an approximately mirror symmetry, here we analyze the mechanism of energy transfer under the duskward IMF as an example. www.nature.com/scientificreports/ www.nature.com/scientificreports/ Figure 1a shows the mechanical energy density flux projected on the YZ plane across the magnetopause as viewed from the sunward direction for B y = 5 nT and P d = 3 nPa. Flowing into the magnetosphere is negative, shown in blue, while the red region indicates the energy flowing out of the magnetosphere. The solid black line represents the boundary where the polar cusp region is located in the plane X = X cusp identified by the pressure gradient and the surface of the magnetopause intersect. Thus the region inside the solid black line represents the dayside magnetopause, and the outer side represents the nightside. Figure 1b shows the distribution of mechanical energy with azimuth angle for the dayside magnetopause (red line), the nightside magnetopause (green line), and the entire magnetopause (blue line), respectively. Figure 1c shows the distribution of the energy integrals of the mechanical energy transfer within the following X-axis intervals: X > X cusp , − 10 R E < X X cusp , − 20 R E < X − 10 R E , − 30 R E < X − 20 R E , and − 40 R E < X − 30 R E . Figure 2 is the same as Fig. 1 but for electromagnetic energy transfer.
It can be seen from Figs. 1a and 2a that the symmetry axes of the magnetopause shape rotate clockwise relative to the meridian plane when only the Y component of the IMF exists. When the dipole tilt is zero, the point of the maximum distance between the magnetopause projected onto the YZ plane and the origin is on the z-axis under the north-south and radial IMFs. However, in the dawn-dusk IMFs, the line between this point and the origin deflects the z-axis by 18 • approximately. For the mechanical energy transfer, on the day side of the magnetopause, one can see a mechanical energy flux inflow near the equator and low latitude, but a small area and magnitude of mechanical energy outflow at the high latitude in Figure 1a. On the nightside, there is a large and weak mechanical energy inflow at the middle and high latitudes in the northern and southern hemispheres, in which a small area of relatively large mechanical energy inflow occurs behind the cusps. At the same time, there is a weak outflow energy flow near the low latitude magnetotail. As shown in Figure 1b, we find that the distribution of mechanical energy with azimuth on the entire magnetopause is not symmetrical between the northern and southern hemispheres but is distorted. It is easy to see that the area with the largest mechanical energy input is near the equator on the day side. In contrast, the mechanical energy input is relatively small at a high latitude. For the night side, the maximum energy input appears near azimuth 75 • and 255 • , i.e., the maximum mechanical energy input occurs on the dusk side in the northern hemisphere and the dawn side in the southern hemisphere, and the mechanical energy output is presented near the low latitude with the maximum output near the equator. In general, the energy transfer on the night side has a relatively large contribution to the mechanical energy transfer across the magnetopause. Figure 1c illustrates that the net input of mechanical energy reaches the maximum on the dayside of the magnetopause, then follows in the near-Earth magnetotail, while the input decreases significantly in the three regions of Figure 2a shows a small area and magnitude of electromagnetic energy outflow at the front of the polar cusp, and the maximum value of the outflow energy flow is around the azimuthal angles 100 • and 280 • . In addition, the rest of the positions on the dayside shows the inflow of electromagnetic energy, and the maximum value of the inflow is located near the equator. Both the inflow and outflow extremes are small and the density flux intensity is weak. However, a large area of electromagnetic energy inflow appears in the far magnetotail with a maximum value near the polar cusp, except in the mid-latitude region of the far magnetotail, where there is a small area and weak electromagnetic energy outflow. Figure 2a shows that the electromagnetic energy flux enters into the magnetosphere in the region near the equator on the dayside, which is confirmed by the distribution www.nature.com/scientificreports/ of the dayside electromagnetic energy integral with an azimuth in Fig. 2b. Meanwhile, as shown in Fig. 2b, the electromagnetic energy input is also present near the azimuthal angles 120 • and 300 • , which is opposite to the outflow shown in Fig. 2a, because there is a small area of outflow but a more extensive area of inflow on the dayside. For the night-side magnetotail, Figure 2b illustrates that the electromagnetic energy input is the greatest in the night-side magnetotail around the azimuths of 135 • and 315 • . In addition, it can be seen that the contribution of energy transfer on the nightside to the entire magnetopause is more significant. Also, the net input of electromagnetic energy is the largest in the near-Earth magnetotail ( −10R E < X X cusp ), which is much larger than the maximum value of the net mechanical energy input as shown in Fig. 2c. Next, we analyze the energy transfer mechanism for duskward IMF. The magnetic field line distribution at B y = 5 nT and P d = 3 nPa is shown in Fig. 3. The gray vector line represents the distribution of magnetic field lines in the YZ plane where X = X cusp is located, the red solid line represents the XZ plane projection of the magnetopause, the black arrows represent the direction of solar wind velocity on the magnetopause surface ( v ), the direction of magnetic field ( B ) is indicated by the green arrows, and the Poynting vector ( S ) is represented by the blue arrows. It can be seen that the reconnection occurs near the dusk side of the northern hemisphere and the dawn side of the southern hemisphere. After the reconnection occurs in both the northern and southern hemispheres, the reconnected magnetic lines are dragged towards the dawnward and duskward at the north and south poles, respectively. At the same time, due to the effect of the solar wind, the reconnected magnetic lines move tailward and eventually become open in the magnetotail. Therefore, the magnetic reconnection occurs on the dusk side of the northern hemisphere and the dawn side of the southern hemisphere, resulting in the accumulation of magnetic field lines near the dawn-north side and the dusk-south side, and the magnetic field lines move toward the magnetotail under the dragging effect of the solar wind.
During duskward IMF, the northern-dusk and the southern-dawn regions of the magnetopause are opened by magnetic reconnection. The accelerated plasma flow is then deflected from the nominal shocked solar wind flow (Fig. 3). As a result, the work of mechanical energy makes the plasma flow enter the magnetosphere from the magnetic field lines opened by magnetic reconnection, along with the magnetic field lines convection to the magnetotail, becoming a part of the plasma mantle, as shown in Fig. 1a. After the asymmetric reconnection of the northern and southern hemispheres, the IMF and the Earth's magnetic field are no longer parallel. According to formulas E = B × v and S = 1 µ 0 E × B , vector B points to the inner side of the magnetosphere, which means that the electromagnetic energy input occurs in the dawn side of the northern hemisphere and the dusk side of the southern hemisphere, as shown in Figs. 2a and 3.
As reported by Lu et al. 2,14 , magnetic reconnection is not the only factor affecting the energy transfer, and viscous effect in the magnetotail also has an important role for the energy transfer across the magnetopause, where the plasma has not only the radial velocity flowing to the magnetotail but also the tangential velocity flowing along the magnetopause from high latitude to low latitude near the magnetopause. To analyze the viscous effect, the velocity vector diagrams (gray vectors) of the streamlines near the magnetopause in the YZ plane for X = −35R E are shown in Fig. 4, where we take the region of Y > 0 and Z > 0 as an example, and the black arrow indicates the direction of the connection between the origin and the magnetopause in the YZ plane, θ is the angle between the black arrow and the positive Y-axis (azimuth angle), and β is the angle between the velocity vector projected onto the YZ plane and the Y-axis. We define β − θ as the streamline deviation angle in the YZ plane. In the absence of external force, the viscous force is zero, and the plasma velocity does not deflect, and then the direction of the black arrow is the direction of the streamline. However, when the viscous force is non-zero, if www.nature.com/scientificreports/ β > θ , as shown in the green line of Fig. 4, the streamline rotates counterclockwise along the magnetopause under the effect of viscosity, otherwise, the streamlines rotate clockwise along the magnetopause, as seen in the red line. The greater the absolute value of the deflection angle, the greater the degree of turning. Figure 4 shows that the deflection direction of the streamlines in the magnetotail is different from that in the north-south and radial solar wind conditions 2,14 . Under the duskward IMF, the streamlines flow from the dusk side of the northern hemisphere (dawn side of the southern hemisphere) to the southern hemisphere dusk side (northern hemisphere dawn side), respectively. The plasma flowing into the lower latitudes accumulates in the southern hemisphere dusk side and northern hemisphere dawn side regions with increased density, in which a small amount flows into the magnetosphere at the lower latitudes, while most of the remaining flow out of the magnetopause. This explains the deflection of the plasma streams at the outflow and inflow positions of the magnetotail under the duskward IMF, as shown in Fig. 1a.

Comparison between dawn-dusk, north-south, and radial IMFs
The magnetopause energy transfer in the the case of radial, the dawn-dusk, and the north-south IMFs are compared. The mechanical energy ( W m ), electromagnetic energy ( W p ), total energy injected into the magnetosphere ( W g ), total energy of solar wind impacting the magnetopause ( W SW ), and energy transmission rate are calculated and compared for the conditions listed in Table 2. As the strength of the magnetic field for increasing duskward IMFs, there is a significant increase in electromagnetic energy and total energy input, while only a tiny increase in mechanical energy. Moreover, the total energy of the solar wind energy impinging on the magnetopause appears to decrease significantly because the surface area of the magnetopause decreases due to the enhanced magnetic field. Figure 5 compares the results of W m , W p , and W g for six simulation cases (radial IMF, dusk-dawn IMFs, and north-south IMFs, respectively), we find that the electromagnetic energy is about twice the mechanical energy input for the dusk-dawn IMFs. For mechanical energy, the inputs are comparable for the radial, dusk-dawn, and north-south IMFs, while for electromagnetic energy, the input in the dusk-dawn IMFs is twice that in the radial IMF, four times that in the northward IMF, but only half that of the southward IMF. The energy transfer rate is then calculated from the total energy W g injected into the magnetosphere and the solar wind energy W SW impinging on the magnetopause, given in the last column in Table 2. The energy transfer rates of 3.58% for dawnward IMF and 3.52% for duskward IMF are higher than 2.22% in the radial IMFs and 1.71% during northward IMF, but lower than 4.95% during southward IMF. Figure 6 shows the energy input as a function of duskward IMF for the simulations with P d = 3nPa . The red line represents the ε parameter, an energy coupling function proposed by Akasofu et al. 43 , and the light blue line represents W g . It can be seen that the total energy injected into the magnetosphere under dawn-dusk IMFs follows the same trend as the energy coupling function proposed by Akasofu et al. 43 , similar to the trend in the north-south condition 2 , but different from the radial condition 14 . However, ε obviously underestimates the energy injected into the magnetosphere. The difference between the formula and the result may be due to the identified magnetopause in this work is an open surface and the plane at the magnetotail is not considered (the application of the numerical simulation results does not allow to consider the energy transfer in the plane of the magnetotail for the time being). The formula represents the energy dissipated in the magnetosphere and ionosphere, which means that not all of it enters the inner magnetosphere, but a portion flows out of the magnetotail. The outflow from the magnetotail is not considered in this work.
In addition, we also compare the deflection degree of the magnetotail streamlines for the dawn-dusk, the radial, and the north-south IMFs (|B|=5nT) under the dynamic pressure is 3nPa, the velocity is 500 km/s, and the density is 7.2 cm −3 , as shown in Fig. 7. Lu et al. 14 have already mentioned in detail that the viscous effect in radial IMF is weaker than in northward IMF. It can also be seen from Fig. 7 that the absolute values of the maximum streamline deflection angle in the magnetotail under the dawn and dusk IMFs are equivalent, and larger than under the radial IMF, but smaller than in the northward IMF.

Summary and conclusions
In summary, the interaction between the solar wind and magnetosphere is investigated by using a three-dimensional adaptive MHD model. The magnetopause surface is identified by the streamline method, and the total energy across the magnetopause under dawn-dusk IMFs is divided into mechanical energy and electromagnetic energy. Magnetic reconnection and viscous interaction processes are used to understand the physical mechanism of energy transfer from the solar wind to the magnetosphere. The main conclusions of this paper are as follows: 1. The symmetry axes of the magnetopause rotate clockwise relative to the meridian plane when only the Y component of the IMF exists. The line between the point on the magnetopause projected on the YZ plane farthest from the origin and the origin deviates from the z-axis, about 18 • under the considered condition in the paper. 2. For the duskward IMF, the energy transfer distribution is no longer symmetric in the north-south direction when viewing from the sunward direction, the mechanical energy flow enters the equatorial and low-latitude regions on the dayside. The electromagnetic energy has a small magnitude of inflow at the equator and low latitudes of the dayside, and a significant energy inflow occurs on the dawn side of the northern hemisphere and the dusk side of the southern hemisphere near the polar cusp. The characteristics of energy distribution across the magnetopause under dawnward and duskward IMF are mirror symmetric about the Y = 0 plane. 3. For the same IMF strength (5nT here), the electromagnetic energy input under the B y oriented IMF is twice as large as the mechanical energy under the IMF B y and the electromagnetic energy under the radial IMF, which is five times as large as the electromagnetic energy during the northward IMF, but only half as large as the electromagnetic energy under the southward IMF. The mechanical energy input under dawn-dusk IMFs is the same order of magnitude as that under radial and north-south IMF. 4. For solar wind condition we considered, the energy transfer rates of 3.58% for dawnward IMF and 3.52% for duskward IMF are higher than 2.22% in the radial IMFs and 1.71% during northward IMF but lower than 4.95% during southward IMF. 5. With increasing IMF, the change of mechanical energy input is insignificant for the dawn-dusk IMFs, while the electromagnetic energy input increases nonlinearly. The Akasofu-type energy-coupling formula, ε , obviously underestimates the energy input from solar wind to magnetosphere for the cases dominated by IMF B y . Figure 5. Comparison of W m , W p , and W g of six simulation cases (sunward, earthward, duskward, dawnward, northward, and southward IMF, respectively). The magnitudes of the magnetic field and dynamic pressure are 5 nT and 3 nPa, respectively.

Data availability
The data comes from the simulation for the steady states under different cases, which depends only on the solar wind conditions for fixed ionospheric conductivities. The datasets used and/or analyzed in this study are available from the corresponding author on request by readers.